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We study the synthesis of lithium isotopes in the hot tori formed around stellar mass black holes 
by accretion of the companion star. We find that sizable amounts of both stable isotopes 6 Li and 
7 Li can be produced, the exact figures varying with the characteristics of the torus and reaching as 
much as 10 -2 Mq for each isotope. This mass output is enough to contaminate the entire Galaxy at 
a level comparable with the original, pre-galactic amount of lithium and to overcome other sources 
such as cosmic-ray spallation or stellar nucleosynthesis. 



INTRODUCTION 

Lithium is one of the very few elements in our Universe 
whose abundance is significantly affected by big bang nu- 
cleosynthesis (BBN), nuclear spallation processes in cos- 
mic rays, and stellar nucleosynthesis. In particular, its 
7 Li isotope is observed in the atmosphere of relatively 
cold stars in our Galaxy, namely halo stars of varying 
metallicity and stars belonging to globular clusters. The 
very first observations pQ have revealed a remarkable flat 
trend in the abundance of 7 Li vs the metallicity of the 
target star, known as the "plateau" ; the most recent ob- 
servations of halo stars (see [5] and references therein) 
show instead a clear "meltdown" of such single abun- 
dance value at the lowest metallicities, as well as a mod- 
erate, yet non-negligible dispersion at intermediate ones. 
The "roof" of this plateau still sits at a factor ~3 be- 
low the prediction from standard BBN, a discrepancy 
that is known as the "lithium problem" . If the value of 
the plateau is of primordial origin, the observations are 
almost impossible to reconcile with standard BBN, and 
exotic models are to be invoked [3 0] . If that value is 
the result of processes of stellar origin, it is possible to 
set up a stellar mechanism able to deplete the original 
abundance [5], but it is still left to demonstration how 
this mechanism can act over three orders of magnitude 
in the metallicity of the star, thus reproducing an almost 
flat envelope. Cosmic-ray spallations can also destroy 
7 Li, but it seems unlikely that they can reproduce the 
observed features [5J. 

The aim of this Letter is not to propose a solution to 
the so-called "lithium problem", but if possible to add 
another piece to the puzzle, as our results go in the di- 
rection of harshening the discrepancy between the BBN 
value and the one observed in stars. This opens up a 
scenario that needs to be taken into account in future 
studies concerning this subject. The existence of stellar 
mass black holes accreting from companion stars is es- 
tablished through observations of so-called microquasars, 
visible in X-ray frequencies in our Galaxy [7]. We show 



that material accreting onto a stellar mass black hole 
from a companion star, in a disk with characteristics of 
a hot torus, can undergo nucleosynthesis, and produce 
sizable quantities of 6 Li and Li. The amount of both 
isotopes depends on parameters such as the mass of the 
black hole, the viscosity of the medium, the accretion rate 
and the expelled fraction, and it can reach up to 10~ 2 M 
within the allowed region of the parameter space. 

We first describe the structure of accretion tori around 
microquasars, and then discuss our implementation of 
the nuclear reaction network and nucleosynthesis equa- 
tions. We finally present our results, as well as draw some 
guidelines and suggestions for future work. 

STRUCTURE OF HOT ACCRETION TORI 

Accretion onto black holes and neutron stars has long 
been recognised as a very efficient means of converting 
rest-mass into radiation. The simplest arguments are 
actually adequate to qualitatively account for the lumi- 
nosities of compact X-ray sources (see e.g. [S] for a re- 
view); yet, self-consistent models are extremely complex 
and have been constructed in only a few limiting cases 
[5Mr7j . Keplerian accretion disks [HI [TO] are amongst the 
most widely-studied models. They form around the black 
hole in a binary as material from the companion star is 
captured and confined to almost-Keplerian orbits given 
its high angular momentum. The viscosity of the orbit- 
ing gas is supposed to play a central role in the accre- 
tion process, but its nature and magnitude are not well- 
understood. Ever since the seminal paper by Shakura & 
Sunyaev [9], it has been customary to parameterise the 
viscous stress in the azimuthal direction as a fraction 
of the gas pressure P, i.e. = aP, where a < 1 is an un- 
known parameter assumed to be constant. The effect of 
viscous stress is two-fold: it removes the angular momen- 
tum at a slow but steady pace leading to a shrinking of 
gas orbits, and it generates (frictional) heat. In standard 
Keplerian disks, electrons and ions in the accreting gas 
are tightly coupled; since the former dissipate energy effi- 
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ciently, the frictional heat induced by viscosity is readily 
radiated away giving rise to a thin and cool disk. 



NUCLEOSYNTHESIS IN HOT ACCRETION 
TORI 



Another type of models, the so-called hot accretion 
tori, was originally proposed in [T3] to account for ra- 
dio jets in some active galactic nuclei. Consider a black 
hole of mass Mbh accreting at a rate M through a 
disk of half-thickness h(r). As in Keplerian disks, the 
momentum transfer is driven by a viscous stress of the 
form /<i = aP, which leads to a radial inflow timescale 
i r (r) = a- l (r/h) 2 (r/rgf^r g /c, where r g = GM BH /c 2 . 
The core idea behind the hot torus model is to conceive 
an accretion configuration in which the inflow timescale 
t r is short enough as to prevent electrons and ions from 
coupling. This is achieved for low accretion rates, specif- 
ically [T3] m = Mc 2 /L Edd < 50a 2 , where L Edd — 
AitGMBH'm'pC/ '<jt is the Eddington luminosity. In such 
regime the frictional heat is not efficiently radiated away 
but is stored, the accretion region is inflated to an ion 
pressure supported torus with h(r) ~ r and ions remain 
at the virial temperature, k E T — (r g /r)m p c 2 /3. It is 
therefore not unnatural to reach ion temperatures of tens 
of MeV in the hot torus. Moreover, the conservation of 
mass combined with the inflow timescale defined above 
gives an estimate of the mean mass density of the gas in a 
cylindrical shell of radius r and half-thickness h ~ r [2] : 
p{r) = a _1 m(r/r 9 ) _3 / 2 3m p /(2(T T r 9 ). Following [T3] we 
take the hot torus to extend from T{ = 100r g down to 
rj = 6r g . Over this range electrons are relativistic [T3] 
and we shall assume they are adiabatically heated inside 
r < r u T e (r) = T(r l )(p(r)/p(r 4 )) 1/3 < T(r). Notice that 
the hot torus is optically thin and not in local thermody- 
namic equilibrium - the gas density is in fact lower and 
the temperature higher than in standard thin disks, sup- 
pressing free-free opacity. The main sources of photons in 
the torus are electron radiating processes. Although syn- 
chrotron and synchrotron self-Compton may play a role 
depending on the magnetic field and electron tempera- 
ture, we shall assume bremsstrahlung is the dominant 
emission mechanism. We make nonetheless an explicit 
test of this hypothesis later on. The photon number den- 
sity per unit energy in a cylindrical shell of radius r, 
dn-y/dE 7 (E 71 r), is computed by taking into account the 
bremsstrahlung emission in the whole torus (for simplic- 
ity we assume spherical symmetry in this step). This 
gives an estimate of the local photon spectrum, which 
will be used in the following. 

It is not straightforward to estimate the fraction of 
microquasars that exhibit a hot torus and such estimate 
would probably span orders of magnitude. However, if 
the accreting material is heated up - or supplied hot - 
and is unable to radiate efficiently, a pressure-supported 
torus will naturally form. Since these requirements are 
astrophysically plausible [13] . we anticipate that a non- 
negligible fraction of microquasars will host hot tori. 



It is clear from the discussion above that extremely 
high temperatures are reached inside the torus and all 
sort of nuclear reactions can take place. In particular, 
reactions with threshold of tens of MeV - unimportant 
in environments such as BBN or the cores of main se- 
quence stars - are easily triggered. This is for instance 
the case of a + a reactions, especially important for the 
synthesis of light elements as pointed out in [T4HT6] . In 
order to pinpoint all the potentially relevant reactions 
for the production or destruction of 6 Li, 7 Li and 7 Be, we 
make use of the network generator NetGen [T7] and the 
NNDC database [TS]. Neglecting processes with at least 
one rare or fastly-decaying reactant or more than two re- 
actants, we identify 74 distinct reactions of which 17 are 
photodissociations. For each reaction the Q-value and 
the threshold are computed with QCalc [TH] while the 
cross-section as a function of the kinetic energy in the 
centre-of-mass frame a(T cm ) is constructed using data 
from NNDC [18] or parameterisations available in the 
literature [2~0t-f2"2"] for specific cases. Whenever available, 
the astrophysical S and R factors are used as well. The 
rate of a given photodissociation is then 



r(r) = j dE 1 ^(E„r) ! j(E 7 )c , (1) 



while for the remaining reactions the thermally averaged 
cross-section at temperature T reads 



(av)(T) = dT cm f(T cm ,T)a(T cm )v(T cm ) , (2) 



with / the Maxwell-Boltzmann distribution, v(T cm ) ~ 
(2T cm / fi) 1 / 2 and (i the reduced mass of the reactants. 
We have cross-checked our (av) (where possible) against 
the parameterised results of [2"lTf2"3"] as well as the pub- 
licly available NACRE database [21] and found good overall 
agreement for the temperatures of interest in those ref- 
erences, i.e. T < 10 9 - 10 10 K. 

We are now in position to write down and solve the 
Boltzmann equation for the abundances of 6 Li, 7 Li and 
7 Be in the hot accretion torus. At every time in our 
calculation we check that the abundances of the pro- 
duced 6 Li, 7 Li and 7 Be are small compared to those of 
the reactants, thus justifying the choice of not following 
self-consistently the abundance of the latter. Consider 
a cylindrical shell of radius r and half-thickness h ~ r, 
shrinking towards the black hole at a pace set by the 
inflow timescale t r . In particular the time-radius conver- 
sion reads t(ri, r) = t r (ri) — t r (r). The mean 7 Li number 
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density in the shell, m Li , obeys 

~~j!r = X! (- 1 + S A l B,y 1 (vv) i n Az n Bi 

i 

where the sums in i, j and k run over processes that 
produce 7 Li, spallation reactions of 7 Li and photodisso- 
ciations, respectively; Aij and Bi are generic reactants 
and SaiBi is unity in the case of identical reactants and 
zero otherwise. Analogous equations hold for 6 Li and 
7 Be. The reactant densities are assumed to track the 
gas density p(r) at fixed abundances (see next Section 
for further details). For a given set (Mbh, M,a) we 
solve the three coupled differential equations between 
(n,ti) = (100r s ,0) and (r f ,t f ) = (6r g ,t(n,rf)) with 
initial conditions ri7 Li = (Yrj J - 1 /mTL l {)p( r i)i where l7 Li 
is the initial mass fraction of 7 Li, and likewise for 6 Li 
and 7 Be. We can then compute the mass of 7 Li in 
the (steady-state) accretion torus at any given moment, 
M* Li = rriTu J r % dr 2nr 2h(r) m^r). After an interval 
of time At, the torus will have accreted a total mass 
M acc = MAt and will have expelled a 7 Li mass 

M t 

Mr Li = f exp -^MAt , (4) 

with Mt the torus mass and f exp the fraction of expelled 
material. Analogous expressions hold for 6 Li and 7 Be. 
We adopt At — 100 Myr and f exp — 0.5 all along, but 
our results can be trivially scaled to other values. A 
value of order unity for the expelled fraction f exp has 
been frequently used [H] [TS] as the presence of powerful 
jets - both observationally confirmed and theoretically 
expected [12] - can lead to significant ejection in hot tori. 

RESULTS AND DISCUSSION 

We have solved the set of coupled Boltzmann equa- 
tions in a large region of the torus parameter space 
(Mbh, M,a). In our default case, we set the initial 
abundances of light elements to the ones leftover by 
BBN, namely gj D/H = 2.49 x 10~ 5 , 3 He/H = 10~ 5 , 
Y P = 0.2486, 7 Li/H = 5.24 x 10~ 10 and 6 Li/H=10~ 14 . 
However, we show later that the initial relative abun- 
dances are not crucial for the final amount of lithium 
isotopes synthesised. The results for sample black hole 
masses of 5 M and 30 M are shown in figure [T] on 
the plane a vs M. It is clear that the nucleosynthesis 
in the hot torus can produce up to > 10~ 4 M Q of 7 Li 
(plus 7 Be) and 6 Li, with little dependence on the black 
hole mass (the small differences between the two cases 
arise mainly because of the different size of the torus and 



inflow time). On the other hand, the Eddington limit 
L = eMc? < LEdd - indicated by the shadowed area in 
figure [l] with typical black hole accretion efficiencies [12] 
0.06 < e < 0.42 - does depend on M BH and sets the 
most stringent constraint on the achievable amount of 
isotopes. We show in figure [2] the total produced mass 
of 7 Li (plus 7 Be) and 6 Li as a function of Mbh for fixed 
a and m. Masses as large as 10 -2 M Q are reached for 
hundred solar masses black holes. For the case of a 5 
M black hole (left panel in figure [T]), the lithium is pro- 
duced at a ratio 6 Li/( 7 Li+ 7 Be) ~ 0.25-0.29. Our results 
agree well with [T5] when we adopt its equation (7) for 
the same black hole accretion configuration; however, in 
figures [l] and [2] we show lithium masses according to our 
equation Q that implicitly assumes a uniform ejection of 
torus material. Computing the expelled mass as if it were 
ejected mostly from the inner regions yields final masses 
higher by a factor < 10, which would only strengthen our 
conclusions. 

The nuclear reactions leading the synthesis of the 
lithium isotopes are the a + a reactions: 4 He(a, n) 7 Bc, 
4 He(a,p) 7 Li, 4 He(a, d) 6 Li and 4 He(a, npfU. All 
present similar thermally averaged cross-sections at the 
temperatures of interest and yield similar outputs of 6 Li, 
7 Li and 7 Be as shown in figure [l] Therefore, the final 
abundance of lithium is sensitive mostly to the content of 
4 He present in the plasma. We find that in order to sup- 
press lithium production by a factor ~10 4 it is necessary 
to suppress the original abundance of 4 He by a factor 10 2 
(notice that the helium abundance enters quadratically 
in equation Q). Varying the initial abundances of D, 3 H 
and 3 Hc negligibly affects the final results, and adding an 
initial abundance of neutrons up to half of the number 
fraction has no significant effect. These checks show that 
the synthesis in the hot torus is sensitive mostly to its 
physical properties rather than the initial conditions of 
the gas chemistry (modulo the abundance of 4 He, which 
is the fundamental reactant for our purposes) . It is worth 
stressing that synthesis at this level is possible because 
of the short resilience time of elements inside the hot 
torus. At photon densities below p/m p , such as the ones 
considered here, the timescale for photodissociation of 
He in all possible channels (that require non-negligible 
amounts of photons above the threshold of ^20 MeV) 
is much larger than the inflow timescale: for the most 
extreme cases we get T~ 1 /t r > 3000. Analogously, the 
photodissociation of 6 Li, 7 Li and 7 Be is also found to be 
unimportant since T~ 1 /t r > 70. Together with the low 
lithium/beryllium densities, this prevents photodissocia- 
tion from being competitive to the a + a reactions. If 
we artificially modify the inflow time, we do in fact ob- 
serve consistent depletion of lithium isotopes for longer 
inflow times, and enhanced final abundances for shorter 
resilience times. Furthermore, we find that scaling the 
photon density by a factor ranging from 10~ 3 up to 10 3 
has little effect on the lithium synthesis, whereas above 



4 




o 





7 Li 






— 7 Ii+ 7 Be 






"Li 






e»i =0.1 






0.01 








or=10-°- 5 






/exp=0.5 






At =100 Myr 



1 10 100 



M BH [Af ra ] 

FIG. 2: The mass of 7 Li (plus 7 Be) and 6 Li synthesised as a 
function of Mbh for fixed a and m. Here we take e = 0.06. 

10 4 there are substantial variations. This large photon 
enhancement is rather unlikely if the photons are to be 
originated in the torus itself; however, if the cool disk sur- 
rounding the torus is nearby and produces large amounts 
of soft photons, these can be upscattered in the torus and 
induce the spallation of light elements as found in [25] . 

The results in figures [l] and [2] are interesting under 
many regards. It is to be reminded that the amount of 
7 Li and 6 Li produced by primordial nucleosynthesis, ac- 
cording to the most recent estimates of baryon density, 
are 7 Li/H = 5.24 x HT 10 and 6 Li/H = 10~ 14 ; this repre- 
sents the minimal "background" level of lithium isotopes 
which are present in the Galaxy at the time of its forma- 
tion. A few simple estimates clarify the relevance of the 
results presented in this Letter: 

• It is expected that 10 8 — 10 9 stellar mass black holes 



are present in our Galaxy, see e.g. [8] and references 
therein. The Milky Way mass in stars and gas is 
Mg a i ~ 5xl0 10 M Q [55], which means that in prin- 
ciple, if a typical hot torus produces 10 -4 M Q of 
both 7 Li and 6 Li, then only 0.1-1% (KT 6 -10~ 5 %) 
of the microquasars in our Galaxy need to host a 
hot torus so that the synthesised amount of 7 Li 
( 6 Li) is comparable to the BBN "background" level. 

• The first generation of stars (Population III, 
Pop. Ill) is expected to form in small halos of to- 
tal baryonic mass 10 5 — 10 6 M Q . Recent numer- 
ical simulations of early star formation show evi- 
dence for multiple Pop. Ill systems, with fragments 
of masses ranging in the sub-hundred solar masses 
regime [27J EH] • It has already been proposed that 
such systems may give rise to microquasars [29] : 
what we wish to highlight here is that the produc- 
tion of 10~ 4 - 10~ 3 M of 7 Li and 10~ 9 - 10~ 8 M 
of 6 Li by a single microquasar is enough to equate 
the BBN "background" level in the halo where the 
Pop. Ill system has formed. These numbers are in- 
deed achieved within the parameter space scanned 
in our study. 

The previous examples are defective in some regards. 
In order for the quoted numbers to be taken at face 
value the material synthesised in the hot torus and ex- 
pelled must be assumed to efficiently mix in the interstel- 
lar medium. This might be the case of primordial star- 
forming haloes - where extremely efficient internal mix- 
ing is expected [3D] -, but not of the Milky Way, which 
raises the issue of even higher concentrations of lithium 
isotopes locally, around the microquasars themselves and 
in regions where microquasars are more frequent. Also, 
although it is not unreasonable that a significant part 
of the material gets expelled from the torus as argued 
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above, presumably some of the expelled material would 
cross cool regions where significant spallation may occur. 
Finally, light element spallation within the torus can be 
important if a nearby, intense external source of photons 
is present. With these caveats in mind, and the need to 
explore further the physics of this "aftermath" part in a 
self-consistent way, it is clear however that the mecha- 
nism described has the potential to produce amounts of 
lithium comparable to the primordial values as well as to 
stellar nucleosynthesis and cosmic-ray synthesis, see dis- 
cussion in e.g. [B]. We point out that future population 
studies should take this mechanism into account, address- 
ing especially the open issues discussed in this Letter. 
It is also worth remarking that microquasars are indeed 
observed in the X-ray bands in our Galaxy, and that 
in principle their optical and infrared counterparts are 
observable. This should give rise to the possibility of ob- 
serving the amount of lithium (whose observable feature 
lies at 670.8 nm), thus checking the validity of the mech- 
anism described so far or, alternatively, permitting to set 
constraints on the parameters and nature of the accretion 
onto stellar mass black holes. Interestingly, anomalously 
high lithium abundances have been observed in late type 
stars, companion of black holes or neutron stars (see [25] 
and references therein), thus suggesting another obser- 
vational strategy to constrain the mechanism we have 
studied. 

We have shown that nucleosynthesis of light elements 
taking place in the hot torus of accreting stellar mass 
black holes can produce sizable amounts of 6 Li and 7 Li 
for a wide range of parameters of the black hole and 
companion star system. Future studies concerning the 
"lithium problem" should take this mechanism into ac- 
count, especially since its actual magnitude and impact 
can in principle be tested through observations. 
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